Faddeev description of two-hole one-particle motion and the single-particle spectral 
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The Faddeev technique is employed to address the problem of describing the influence of both 
particle-particle and particle-hole phonons on the single-particle self-energy. The scope of the few- 
body Faddeev equations is extended to describe the motion of two-hole one-particle (two-particle 
one-hole) excitations. This formalism allows to sum both particle-particle and particle-hole phonons, 
obtained separately in the Random Phase Approximation. The appearance of spurious solutions for 
the present application of the Faddeev method is related to the inclusion of a consistent set of dia- 
grams. The formalism presented here appears practical for finite nuclei and achieves a simultaneous 
inclusion of particle-particle and particle-hole phonons to all orders while the spurious solutions are 
properly eliminated. 

PACS numbers; 21.10.Jx, 21.60.-n, 21.60. Jz 



I. INTRODUCTION 



In recent years, the study of (e, e'p) reactions has been one of the most useful tools to probe correlations in nuclei. 
Absolute spectroscopic factors have become available for many closed-shell nuclei ||l|-^ demonstrating that the removal 
^ . probability for nucleons from these systems is reduced by about 35% in comparison with the simple shell model. The 
theoretical description of this reduction requires the inclusion of both short-range and long-range correlations. For 
nuclear matter a strength removal of about 15% is obtained by including short-range correlations For ^^O the 
inclusion of short-range correlations leads to removal of single-particle (sp) strength of the order of 10% The 
T-H inclusion of long-range correlations for heavier nuclei like ^*Ca yields a qualitative description of the sp strength 
distribution by including in the nucieon self-energy the coupling to either low-lyingcoUective particle-hole (ph) or 
' particle-particle (pp) phonons calculated in Random Phase approximation (RPA) The additional depletion of 
(~| I about 10% due to short-range correlations for this nucleus leads to a reasonable quantitative agreement for the largest 
' fragments of the experimental strength distribution. 
'T^ . Theoretical calculations of hole spectroscopic factors for ^^O are not so successful. The experimental spectroscopic 
^ ' strength Q for the knockout of a proton from both the a-nd P3/2 shells corresponds to about 60%. The 10% 
(— I reduction due to short-range correlations is mostly compensated by the proper inclusion of the center-of-mass motion 
j^' ■ which enhances the probability for p removal by about 7% |^ . Calculations based on the Green's function approach, 
including both long- and short-range correlations, yield about a 25% reduction |l^. These results still need to be 
corrected for the center-of-mass effect. It is therefore fair to conclude that the present theoretical results for ^^O are 
' still about 20% away from the experimental data. The importance of low-energy correlations is clearly demonstrated 
by the results of Ref. and their proper inclusion is therefore crucial for a complete understanding of this puzzle. In 
the latter work the self-energy was obtained including the effects of interactions between both pp and ph excitations 
in the Tamm-Dancoff approximation (TDA). In order to account for the coupling to collective excitations that are 
actually observed in ^^O it is necessary to at least consider an RPA description of the isoscalar negative parity 
states |ll[] . To account for the low- lying isoscalar positive parity states an even more complicated treatment will be 
required. Sizable collective effects are also present in the pp and hole-hole (hh) excitations involving tensor correlations 
for isoscalar and pair correlations for isovector states. Another argument to improve the description of the coupling of 
sp states to low-lying collective excitations is provided by the lack of fragmentation at low energy obtained in present 
theoretical studies |10[| in disagreement with experimental data. 

One of the goals of the present work is to account for the collectivity in the ph and pp (hh) channels in a consistent 
way while including these excitations at least at the RPA level. Since the observed fragmentation and depletion of 
the sp strength in ^^O is quite substantial, it is reasonable to assume that these features are also important in the 
description of the excitations that contribute to the self-energy. This results in a self-consistent formulation where 
the dressing of the nucleons is incorporated in the description of the collective excitations that ultimately leads to the 
dressing itself. For this reason the present work will be formulated using self-consistent Green's functions (SCGF). 
This type of self-consistency must also be considered in describing pairing correlations in semi- magic nuclei ||l^,|l^ . 
A formalism in which both pp and ph phonons are treated at the RPA level in the self-energy was proposed 
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in ll^. This work focused on the two-time two-particle one-hole (2plh) propagator and generated a formulation that 
reduces to either including the pp-RPA or the ph-RPA phonons in the self-energy when ph or pp vertices are omitted, 
respectively. This expansion, however, was obtained using some drastic assumptions and disregarding some of the 
constraints that arise when propagators in different diagrams of the expansion are connected. Related to this issue is 
the appearance of unphysical solutions for the 2plh propagator that have been discussed in |15|. Results in have 
therefore been obtained mostly for the TDA treatment of the 2plh propagator. Moreover, this approximation was 
obtained by employing mean-field (single-pole) sp propagators. In the present self-consistent treatment, which sums 
fully dressed propagators, this approach is no longer possible. 

To proceed with the inclusion of both pp and ph collectivity in the nucleon self-energy it is important to note 
that the naive summation of diagrams containing both pp and ph phonons leads to serious inconsistencies. This 
approximation is depicted in Fig. |l|. The last of the three diagrams on the right hand side is already contained in each 
of the other two and must therefore be subtracted to avoid double counting. This subtraction introduces spurious 
poles in the Lehmann representation of the self-energy and generates meaningless solutions of the Dyson equation. The 
minus sign in front this term may also prevent in some cases the proper normalization of the spectroscopic amplitudes. 
This feature can be understood by considering a possible solution near such a spurious pole. The normalization is 
determined by the derivative of the self-energy at this energy JTst and will not yield a correct result on account of the 
additional minus sign when the third diagram dominates. In addition, each of the first two terms in Fig. |l| ignores the 
Pauli correlations between the freely propagating line and the quasi-particles forming the phonons, as noted in In 
the present work a formalism is pursued which sums the contribution of the pp and ph phonons to the self-energy to 
all orders avoiding the subtraction of the second order diagram. The treatment of Pauli correlations is improved over 
methods that employ ph RPA phonons in the self-energy since all exchange terms at the 2plh level are consistently 
included. 

Other approaches have been proposed in the literature that attempt to extend the nature of the phonon correlations 
included by performing massive summations of diagrams p6|"p^. Nevertheless, a consistent resummation of both pp 
and ph phonons to all orders has not been achieved in these papers. The main problem in pursuing such an infinite 
summation of diagrams for the 2plh propagator, which includes both pp and ph RPA correlations, is related to 
the fact that a two-body interaction can invert the sense of propagation of only two lines (i.e. change at most two 
holes in two particles and vice versa) while the third line continues to propagate in the original direction. In this 
way, a propagator depending on more than two times is generated. It is therefore necessary first to consider an 
exact formulation involving the four-time Green's function for the 2plh propagator. Direct application of four-time 
propagator equations presumably will remain impractical for the forseeable future. Appropriate approximations to 
this equation are therefore necessary to construct the relevant two-times Green's functions which contains the sought 
after correlations. The scheme studied in this paper consists in computing the RPA phonons in the pp and ph 
channels, separately, and then summing them to all orders employing a Faddeev technique ||l9| , ^0| . 

A nontrivial problem in the implementation of the Faddeev equations is the appearance of spurious solutions ^1 22 
which also have to be considered for the 2plh propagator. As pointed out in |2^], the spurious eigenstates are easy 
to recognize for the few-body problem since they also diagonalize the unperturbed hamiltonian. Their main features 
are that their eigenvalues are known and that their wave-function amplitudes sum up to zero. The situation is more 
complicated in the many-body problem when the Faddeev technique is employed. In particular, the fulfillment of 
closure relations for pp and ph amplitudes is related to the behavior of the spurious Faddeev eigenstates. Without 
a consistent treatment of this relation the spurious solutions will mix with the physically meaningful ones. Applying 
the Faddeev technique to the many-body problem, it is important to solve for all physical solutions that contribute 
to the self-energy. Thus, it is necessary to develop a formalism in which the spurious solutions are correctly separated 
from the physical ones. 

The practical implementation of the present Faddeev scheme is beyond the scope of this paper. The resulting set of 
equations require a great deal of computational effort, especially when dressed propagators are employed. Nevertheless, 
it appears that they can be solved in practice and results using this formalism will be presented elsewhere. 

In Sec. H we briefly describe the SCGF approach based on the Dyson equation and present the exact Faddeev 
formalism for the four-times 2plh propagator. The construction of a consistent formulation for the two-time 2plh 
propagator including the propagation of the pp and ph RPA phonons to all orders is presented in Sec. III. Although 
sofar only the 2plh propagator has been mentioned, it should be understood that the corresponding two-hole onc- 
particle (2hlp) propagator must be included in the calculation of the nucleon self-energy. In the present work no 
coupling terms are considered that transform the 2plh into the 2hlp propagators (or vice versa). For this reason the 
same technique can be used for both propagators and we will use the generic 2plh to represent both. In Sec. IV the 
appearance and treatment of spurious solutions is discussed. Some technical details are relegated to the appendices. 
Conclusions are drawn in Sec. M. 
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II. SELF-CONSISTENT GREEN'S FUNCTION APPROACH AND 2P1H PROPAGATOR 



A. Self-energy and 2plh propagator 

We consider a finite system of A fermions interacting by means of a two-body interaction V. As usual one may 
introduce an appropriate mean-field potential U to localize the nucleons and split the hamiltonian into an unperturbed 
one-body part Hq = T +U and a residual interaction Hi = V — U . Since we are mainly interested in low-lying bound 
states of finite systems, we consider sp states with discrete quantum numbers. As a basis, we choose the set of sp 
states {a} that diagonalize Hq with corresponding eigenvalues £° . The total hamiltonian can then be written as 

H ^ Ho + Hi ^ cica + i ^<^0>iS 4. cj, cs - '^Ua,p cicp , (1) 

a \ afi'^S af3 j 

where c]j (cq) are the creation (destruction) operators of a particle in the state a, Ki/j.^i are the antisymmetrized 
matrix elements of Y ^ and Vafi correspond to the matrix elements of \J . 

The one-body propagator of the A-body system with ground state l^*;^) is defined as |^ , ^ 

S.Mr) = -*«|r[c„(r)ct(0)] 1*^), (2) 

where cj,(t) and Caii) now correspond to operators in the Heisenberg picture. In the Lehmann representation, all the 
eigenvalues of the excited states of the systems with A -I- 1 and A — 1 particles appear, as well as their spectroscopic 
amplitudes for transitions to those states that are relevant for comparison with experimental data. 

The propagator 5q/3(t) can be obtained as the sum of an infinite set of diagrams, built from interaction vertices 

and unperturbed sp propagators corresponding to Hq- In the nuclear case, a strong coupling exists between the 
sp degree of freedom and both collective low-lying states as well as high-lying states. The latter coupling is related 
to the strong short-range repulsion in the nuclear force. The resulting fragmentation of the sp strength, as observed 
in experimental data, suggests that this feature must already be included in the description of these couplings. For 
this reason self-consistent one-body propagators need to be considered in the construction of the nucleon self-energy. 
This self-consistency feature also emerges in an exact formulation, involving the coupling to two-, three- and A-body 
propagators, which can be derived using the equation of motion method | |2^ . In short, this means that for the 
nuclear case one needs to develop the perturbation theory in terms of the dressed propagator (^) approximated in an 
appropriate way. 

The approach we use here consists in computing gap (t) as a solution of the Dyson equation 

9Mr) = gi°j(r) + .g(") (r - - ^2 ^(is) , (3) 

where S*^(t) is the irreducible self-energy. Here and in the following, we employ the convention of summing over all 
repeated indices and integrate from — oo to -|-c» over all repeated time variables, unless specified otherwise. 

By considering the equation of motion for gapiT), one obtains that S]*^(r) can be written as the sum of two terms 

I]*^(t) = SfJ + V^x,^., R^,x,Mt-,t,t+;0+, 0,0-) V^s,i3e , (4) 

where Ti^^ represents the (time independent) Hartree-Fock part of the self-energy, which can be computed from the 
solution 5q/3(t) itself. The 2plh propagator R, appearing in the last term of Eq. (^), contains the sum of all so-called 
one-particle irreducible diagrams which cannot be separated by cutting a single line. These terms are included in the 
2plh Green's function, g'^Jx^ff-y' defined below. The relation between R and g'^P^'^ is given by |2^ 

Ra (^1,^2,^3! ^4, t^^ta) 

= 3'^.W7(*i>^2,i3;t4,t5,<6) - 9iiM^hM-M,t') g-^{t' ~t'') gl^i^„p{UX-MM) , (5) 

in which g^^ is the inverse of the one-body Green's function (H) and g^^ and g'^P^'^ are the 4- and 6- point Green's 
functions defined as 

g'Jp^^sitiMMM) = -I «| T[cMi2)ca(ii)4(t3)4(t4)] lO (6) 

and 
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4^W7(*l'*2,*3;i4,t5,t6) = -i (*^| T[4(i3)c.(t2K(il)cl(i4)4(i5K(i6)] l^-^) , (7) 

respectively. 

The propagator R^vX.ap-t is the solution of the following equation which has a similar form as the Bethe-Salpeter 
equation for pp and ph propagators 

Ra 

— 9tia{ti ~ t4)gi,i3{t2 — t5)g^x{tQ — ta) — g^aih — t4)gfii3{ti — t^)g^\{tQ ~ t^) 
+gfif,'{ti - t[)g^^>{t2 - t'2)g\'x{t'^ - h) 

X ii'^',y'A',a'/3'7'(^l: ^2: ^4: ''e) -^a'/5'7' ,Q/37 (^41 ^5 1 ^6 i ^4 , ^5 i '^d) ; (8) 

which is shown in Fig. ^in terms of Feynman diagrams. The interaction vertex, also shown in Fig. ^, is given by 

{tl,t2, t4, is, tg) 

= K^Jt:Uit2,h; t5,te)g~^ih - U) + K^^x%iti,h; U, U)g-^p{t2 - h) 

+ U, h)9-l{h - h) + K^Pxi(3^{hMM: U, hM) ■ (9) 

In Eq. (H), K'^pp'^ and K^p^'^^ represent the pp and ph irreducible vertices while K^p^^^ is the 2plh irreducible vertex. 
It should be noted that in Eq. (|^) the propagator R^uX^aPf is only required at two times and therefore its complete 
knowledge, as given by is not necessary to solve the Dyson equation. On the other hand, the dependence on the 
time variables ti, t2 and ^3 is employed in the Bethe-Salpeter Eq. (^), thus requiring that at least a 4-times object be 
employed to solve for the 2plh motion exactly. 

Equations (^), (0) and together form a set of coupled equations, where the same propagator, which solves the 
Dyson equation (g), appears as input in the Bethe-Salpeter equation (p|). If the irreducible vertices K^pp\ k^p'^'> 
and i^(PP^) are also expressed in terms of the (/a/3(T), then Eqs. (§) and ^) will generate a self-consistent expansion. 
Obviously, Eq. (||) and the irreducible vertex (^) represent the exact solution for R and therefore require a suitable 
approximation. 



B. Faddeev-Bethe-Salpeter equations 

Eq. (^) can be reduced to a set of coupled equations in a way similar to the method proposed by Faddeev to solve 
the three-body problem [p^ , ^ . The inclusion of pp and ph RPA phonons in a consistent way requires this Faddeev 
approach since it provides a natural framework for correctly iterating quantities that have already been summed to all 
orders like these RPA phonons. In the present work we will neglect the contribution of the irreducible K'^pp^^ term in 
Eq. since it leads to the coupling of higher order particle-hole terms than already considered in the following. We 
will therefore require only three Faddeev components. Following standard notation in the literature |2^, R'^^Ix apj 
will represent the component related to all diagrams ending with a vertex between legs j and k with {i,j,k) cyclic 
permutations of (1,2,3). We will employ the convention in which the third leg propagates in the opposite direction 
with respect to the first two. The Faddeev components i?'*-* can be written in terms of the 2plh propagator R and 
the contribution of the three dressed but noninteracting sp propagators. This definition is given in detail here for all 
three components, omitting explicit reference to the time variables for convenience of notation 

^^I^^X,al3'y = 9i^e9pX K[p^^^ Rfirja,al3-t + ^ i9pa 9i^0 9lX - gua 9tiP 9ix) , (10a) 
^'puX,ap-y = 9tJ.t9pX K'fp^ Rriucr,aP'i + ^ ^9 pa 9vl3 9^^ ~ 9va 9pp 9jx) , (10b) 
RpuX,al3'y = 9p<i9yp K^fp% RriaX.apj + ^ {9pa 9u0 9lX " 9i^a 9pP 9jx) ■ (lOc) 

The factor i in Eqs. ( p^ ) properly takes into account the exchange symmetry between the parallel lines in the Faddeev 
equations. With these definitions the full propagator (||) is given by 

RpiyX,al3f = ^ Rplx,al3'y ~ pa 9v0 9lX - 9va 9 pP g-yx) ■ (H) 

i=l,2,3 
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The Faddeev equations now take the following form 

^tl^ap-y = 2 ~ ^^'^ ^^^^ 

+ ffMM' ff-' 5A'A T%,^,^^„^„^, {R^X"X".o.f}^ + 4'l"A",a/3^) , * = 1, 2, 3 (12) 

where the ^ ^1/37 vertices obey the following symmetry relations and are defined by 

r^tlA,a/37(^l'*2,t3;i4,i5,i6) = r[,^jA,/3a7 (^2 , tl , ^3 ; ^5 , ^4, ^6 ) = 5/7^(^1 " f (^2 , ^3 ; ^5 , ^e) , (13a) 

r^tA,a/37(*l'*2,i3;i4,i5,i6) = T^^^^^^^^^ (i2 , , ^3 ; ^5 , ^4 , ^e) ^ Q^xih - h) f ^'L^^/J (^1 ' *2 ; ^4 , ^s) ■ (13b) 

The gamma matrices F'^^^') and r(p'') are the four-point functions that solve the Bethe-Salpeter equation for the pp 
and ph motion. These vertex functions contain the pp and ph phonons and can be written as 

^7X0/3(^1' ^2; is, <4) = K'^sl,p{ti,t2]h,'t'i) 

+ ^2; g,,{t[ - t',) g,At'2 - Q K%%(t'^,t'^- h, U) , (14a) 

r75^1/3(*l'*2;i3,i4) = K''^i^lp{ti,t2]t'i,tA) 

+ fS,(ii,i2;i'i,ti) gM-t'^) g.At'i-Q K^^.%(t':iA-MM) ■ (I4b) 



III. APPROXIMATE FADDEEV EQUATIONS FOR 2P1H MOTION 



Apart from neglecting the K^'^p^^'' vertex, Eq. ( |l^ is otherwise a complete equation for the 2plh propagator. 
This general equation involves quantities which depend on several times and is therefore too complex to be solved 
numerically. In order to construct a manageable approximation scheme that includes the relevant physical ingredients 
two simplifications will be considered in this section. The first one involves the restriction to two-time pp and ph 
vertices that include the respective RPA contributions in these channels. This approximation is the minimum step 
that maintains the simultaneous inclusion of both pp and ph collective low-lying excitations in describing the sp 
propagator. Second, it is necessary to simplify Eq. (|l2|) to include only two-times Green's functions. This procedure 
no longer allows the inversion of the the propagation direction of all three lines together. As a result, the Faddeev 
equations split up in two separate expansions for the 2plh and the 2hlp components. Although the hole spectral 
function is of primary interest for comparison with experimental data, it must be stressed that both 2plh and 2hlp 
components are needed to generate the self-consistent solution for the sp propagator. Since the formalism involved is 
the same for both components, we will describe only the forward-going (2plh) expansion. The equations for the 2hlp 
case are completely analogous. 



A. Faddeev equations 



To construct the present approximation scheme, it is more convenient to use the energy representation. The 
corresponding Lehmann representation of the sp propagator (||) is given by 

gaf3{uj) = > + > _ . , (15) 

^ uj-eZ +tr] ^ uj-e,. -17] 

where A"" = i^n'^^l'^^l'^o) (^q — ~^|cq|^o )) spectroscopic amplitudes for the excited states of a system 

with A + 1 (A — 1) particles and the poles e+ = — Eq (e^ = Eq — E^~^) correspond to the excitation energies 

with respect to the A-body ground state. In Eq. ( p^ ) and in the following, we use the indices n and k to enumerate 
the fragments associated with the one-particle and one-hole excitations, respectively. 
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Employing the bare interaction Vap^-yS for the vertices K^pp^ and K^p^\ the Bethe-Salpeter equations (14) reduce 
to the usual dressed RPA (DRPA) equations The solutions of these equations depend only on two times. 

These pp and ph phonons correspond to the dressed version of the phonons that are considered in Ref. (see also 
Fig. [^). These excitations describe the correlations that we aim to iterate to all orders and, subsequently, to include 
in the self-energy as explained in the introduction. These DRPA solutions can then be substituted in Eqs. (13) to 
generate the F*^') matrices to be used in the Faddeev expansion. Both the forward- and backward-going components 
of the DRPA solutions arc included into the expansion as illustrated in Fig. | . This is crucial in order to eliminate 
the spurious solutions of the Faddeev equations as will be explained in Sec. |V and Appendix |^. 

The working expression for the F^'^ matrices, which depends on only two times (or equivalently one energy), is given 
in some detail in Appendix Here we only need to stress that the resulting F^'^'s cannot invert the freely propagating 
line from hole to particle or vice versa, i.e. they cannot connect the 2plh amplitudes with the 2hlp ones. For this 
reason, the pp and ph phonons will be summed only in one time direction in a TDA way contributing separately 
to the 2plh and 2hlp propagators. The reader may notice that two contributions of the type shown in Fig. ^ can 
connect the 2plh and the 2hlp propagators. The inclusion of such terms leads to the simultaneous propagation of two 
phonons which requires an extension of the approximation presented in this paper. Since these terms are expected to 
contribute only in higher order, we will neglect them in the following. We note that the collective RPA correlations 
in the pp and ph channels have already been computed through Eqs. ([ij) and therefore remain properly included in 
our approximation. 

The remaining complication, related to the use of dressed propagators, concerns the interactions vertices (p^). As 
explained in Appendix ^ the F*^*^ and the propagators i?*^*^ need to be redefined in such a way that their matrix 
elements also depend on the indices (n, n', k), which label the fragments of the propagators. This implies that the 
eigenvalue equations will involve summations on both the sp indices (a, /3, 7) and the ones corresponding to the 
fragmentation, («„, n/3, A;^). The 2plh propagator and its Faddeev components, as defined in Eqs. (||) and (|To|), are 
recovered only at the end by summing the solutions over all values of {ua, np, k^) and (n^, Uj/, kx). 

Putting together all the above considerations, the resulting approximation to the Faddeev equations (|l^) can be 
rewritten in a way where all the propagators involved depend only on one energy variable (or two time variables). 
The forward-going part of this expansion can be written as follows 

n>(j) I \ 

2 \^ ^n^^i^n,yXkx,anaf3nfjjk.y\^ J ^ uni,^n^Xkx,anal5ni^^k~f\ } J 
+ ^ vn^^in^Xkx.fi'n'^v'n'^X'k'^i^) ^ u' n^fi' n'^^X' k'^,p."7i'^u"n'J X" k'^(^) 

X fi?^^2 „ „ „ , (W) + R''':} „ „ n.ni.,, a i. (t^)! , ^=1,2,3. (16) 

In Eq. (|l^, G^^ is the forward- going part of the 2plh propagator for three dressed but noninteracting lines. Using 
the notations introduced after Eq. dla) we have 



G u.n^vn^Xkx,an^3n„-ikA^^ — ^n^^.n^ ^n^,np ^kx,kj TTT ; T - n , . ■ (17) 



Eqs. (p^, together with the F^'^'s given in Appendix^ approximate the general "Faddeev-Bethe-Salpeter" expan- 
sion to a tractable set of equations involving only two-time objects. It is important to note that these equations are 
still expressed in terms of the self-consistent solution gap{Lu) and include both pp and ph RPA phonons in a correct 
way. Thus they maintain all the features relevant for the physics we aim to describe. 



B. Faddeev amplitudes 

Eqs. (p^ ) involve the use of propagators depending on a large number of indices. As a consequence, the dimension of 
the problem could easily grow up to a point where no practical application is feasible for a real system. This difficulty 
can be overcome by introducing a new set of spectroscopic amplitudes which depend only on the indices labelling the 
particle and hole fragments (n, n' , k) ||3l|,P4 . Thereby the problem is reexpressed by changing from the basis of sp 
states {a}, used in definitions (^ and (|10|)7 to a new formulation constructed in terms of the fragments labeled by 
{n, k}. This procedure also allows to rewrite the eigenvalue and normalization conditions corresponding to Eqs. (I^) 
in a more concise way. As long as the interaction elements Va/s^-yS are energy independent, all the solutions can then 
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be obtained through a single diagonahzation. This approach is particularly satisfactory from a physical point of view 
since the equations reflect the mixing of the 2plh states represented by the (n, n' , k) fragments. This new formulation 
does not introduce any further approximation. Nevertheless, since it appears relevant for a practical solution of the 
problem, we will describe it in the following. 

The Lehmann representation of the Faddeev components i?^*^ ^n^xkx an^/Snp-yk > contains all the poles e^f of the 
2plh propagator, each with its own residue. One obtains for these components 

m ™ ' 

where the superscript m labels the solutions of Faddeev equations. In Eq. (p^, R^^l^^ represent components containing 

the same poles as (0). The sum of these terms cancels exactly the contribution of the three freely propagating 
lines in Eq. (pi]), leaving in the Lehmann representation of the Faddeev propagator only those poles e^, that 
correspond to correlated 2plh states. This is most easily demonstrated by applying the the DRPA equations to both 
sides of Eq. (llF 



The vectors represent the amplitudes of the three Faddeev components and sum up to the residues of 



the full propagator 



i=l,2,3 



We now define new Faddeev amplitudes x'nin'-ik which are related to the (5^^^ 's in such a way that 

t-'anif}n2-lk ~ ■J'-y -^nm^k ' V^"7 

where no summation is performed over the particle and hole indices ni, n2 and k. We also introduce the notation for 
the spectroscopic amplitude, analogous to Eq. (H) 



1=1,2,3 



In general, Xmn^k and the components x'^^_^^^j, define four vectors x and x^'^ all belonging to the same linear space. 
It is useful to split up the latter in two spaces Va and Vs containing all the vectors that are antisymmetric and 
symmetric with respect to the exchange of the two particle indices ni and 712, respectively. Thus, 



We also define a vector X containing all the three components 



(22) 



X = x(2) ^ VX®VI . (23) 



Here and in the following, we use the convention to denote vectors with lower case boldface and operators (matrices) 
with plain capital letters belonging to the space Va^Vs- Both vectors and matrices in the space ® Vs are denoted 
by capital boldface letters. We will also use I for the identity matrix in the Va ® Vs space and the superscript ex to 
indicate the vectors obtained by exchanging the two particle indices (thus, /^^ is the operator exchanging ni and ?i2 
in Eqs. (|0l) and (|l|) ). 




C. Faddeev Hamiltonian 



The eigenvalue equation for the Faddeev expansion can be obtained by substituting the Lehmann representation ( p^ ) 
into Eqs. ( |l^ ) and extracting the residues of the poles. After some algebra, one obtains the following set of 
equations in terms of the x^'^ vectors 
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rii) 



1,2,3 



(24) 



In Eq. (p4|), the components of the matrices i?^'^ J7''^ and T^'^ are related to the spectroscopic amphtudes of the 
DRPA propagators, as explained in Appendix |b[ The D*^*^ 's are diagonal matrices containing the eigenvalues of the 
corresponding DRPA. 

One can now define block-diagonal matrices H, D, etc... which contain on the diagonal the matrices 77'*^ D^-^\ etc. 
These matrices act on the vectors X defined in (p3|). Using this notation one combine Eqs. (p4|) as follows 



HHt 



U- 



where we have also introduced the matrix 



M 



w D 



/ / 
/ / 



-Tt 



M X 



(25) 



(26) 



that takes into account the proper mixing between the Faddeev components. 
By introducing the vector 



Y = 



1 



-Tt M X 



(27) 



(which appears in Eq. ( |25| ) ) and remembering that D is a diagonal matrix, it is possible to manipulate Eq. ( p5| ) into 
the usual form of an eigenvalue equation 



w X = F X 

where we have introduced the Faddeev hamiltonian F ||2^ , which is given by 

F = [I-HH^M]"^ U {TM + D(U-^) [l-HHtM]} 



(28) 



(29) 



The form (psj ) of the Faddeev eigenvalue equations is useful, since it reduces the problem to the diagonalization of a 
single (non-hermitian) hamiltonian. 

The hamiltonian F can still correspond to a large matrix requiring a large amount of CPU time to diagonalize. As 
will be explained later in Sec. [V, about 2/3 of its solutions are trivial and without physical meaning. Thus, it is not 
necessary to diagonalize the full Faddeev hamiltonian (E9|) but one can project it onto the space of physical solutions. 



D. Symmetry Requirements and Normalization Conditions 



As a consequence of the Pauli exclusion principle, the spectroscopic amplitudes for the 2plh motion have to be 
antisymmetric with respect to the exchange of the two particle indices. This statement applies to the full spectroscopic 
amplitudes ( |l9|) and (|2^) but not to the single Faddeev components, which have more complicated exchange properties. 
To exhibit the correct symmetry requirements for the Faddeev components, it is useful to introduce the following 
exchange operator, which works on the space (ESh of the three x^*^ components: 



(30) 



The form of matrix ( ^0|) takes into account that the component x^^' has to change into x*-^^ when the first two legs 
(i.e. the two particles) are exchanged. Since P is idempotent (i.e. = P), it has only eigenvalues +1 and —1 and 
the respective eigenvectors are of the form 



X_i = 



Tex 

1 J^, 

jex 



Xfc 



and X+ 1 = 



Xa 

+/'^^ x, 

jex 



Xb 



(31) 
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in which and Xb are any two vectors. One easily recognizes that the three Faddeev components of X_i and X+i 
give rise to antisymmetric and symmetric spectroscopic amphtudes, respectively, when they are inserted in Eq. (|2l|). 

Using the symmetry properties of the interaction boxes ( [l^ ) and the definition of M (|2^), one can show that P 
commutes with the matrix multiplying X in Eq. ( ^5|) and therefore with the Faddeev hamiltonian (|2^). Thus, P and 
F must have a common set of eigenvalues. The relevant eigenvectors in the present case correspond to those involving 



R 



('0 



The normalization condition is derived as usual by considering the Lehmann representation (|18|) for the components 
(lu) p3|,p3. One can expand around a given pole and consider terms to order zero and then make use 



of the conjugate of the eigenvalue equation (p4[). The result is a condition for the X^'^'s which only allows proper 
normalization for the antisymmetric component. These antisymmetric solutions X_i satisfy the following condition 



x^x 



E 

1=1,2,3 



(32) 



where x is the spectroscopic amplitude appearing in Eq. ( pT] ) and the factor of 2 appears because a sum over all indices 



of: 



'nin2k 



j) is implied, which includes also the exchange terms. Eq. (|32| ) differs from the usual normalization of a 



wavefunction for the fact that we have to subtract the additional terms 



D(0 



T 



(i)t 







■(i) 



1,2,3. 



(33) 



These contributions correspond to the diagrams shown in Fig. ^ which have been discarded in the present expansion. 



IV. TREATMENT OF SPURIOUS SOLUTIONS 



The Faddeev formalism is based on the introduction of different components x^'^ , which belong to the same linear 
space of the total spectroscopic amplitude x (^2|). These components are the solutions of the Faddeev-eigenvalue 
equation (|28|), which is formulated in a larger space in terms of the vectors X containing all three x^'^. Only one third 
of the solutions in this larger space have physical meaning while the others have to be discarded. One can clarify this 
problem by looking at how the complete spectroscopic amplitudes x are obtained from the components x^*) through 
Eq. (^). Relevant details for treating this issue are discussed below. 

The antisymmetric solutions of the Faddeev equations X_.i are determined from two independent vectors x^ and 
xj as shown in Eq. (pT]). In particular, one has to specify both the symmetric and antisymmetric parts of the first 
vector (xa) and only the antisymmetric part of the seocnd (xt,). These solutions therefore belong to the space 



Va ® Va ® Vs 



(34) 



The complete spectroscopic amplitudes x must also be antisymmetric under the exchange of the two particle indices, 
so they belong to Va- Thus, Eq. ( pT| ) must be regarded as a projection from Vf to the smaller space Va and therefore 
must have a nonvanishing kernel. We denote this kernel by Vsp and refer to its vectors as spurious states^ '^Sp- 
Although these states satisfy the Pauli requirements, they don't yield any contribution to the full 2plh propagator 
and have no physical meaning. We also consider the space Vph, which is orthogonal to the kernel Vsp and contains 
the antisymmetric states Yp^ which generate nonvanishing spectroscopic amplitudes x. The vectors belonging to Vph 
produce contributions to the 2plh propagator and therefore in the following they will be referred as physical states. 
In Appendix |^, explicit basis sets for the Vph and Vsp spaces are given. Obviously the combination of these two basis 
sets forms an orthogonal basis of Vp and one has Vph = Va and Vsp = Va <Xi Ks. 

It must be stressed that in general the solutions of the Faddeev eigenvalue equation (^8|) do not automatically 
separate into the physical and spurious states just defined. Nevertheless, it is shown in Appendix ^ that the states of 
Vsp are proper eigenstates of the Faddeev hamiltonian for the expansion presented in this paper. This feature always 
occurs for the three-body problem but is not guaranteed when working with quasiparticle and quasihole excitations 
unless a proper set of diagrams is considered. When this condition is satisfied, there exist a set of spurious solutions 
of the Faddeev equation ( |28| ) which spans the space Vsp completely. The projection of the Faddeev hamiltonian ( [29| ) 
onto the physical and spurious subspaces Vph and Vsp then takes the form 



(Ph\F\Ph) 
{Sp\F\Ph) (Sp\F\Sp) 



(35) 
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It should be noted that the physical states Yp/j, belonging to Vph, differ from the spurious ones Ysp (€ Vsp) not 
only because they give rise to physically meaningful spectroscopic amplitudes but also because they are not solutions 
of the Faddeev equations (p8|). In general, a physically meaningful eigenvector of 'X.physicah is a mixture of states 
belonging to both Vph and Vsp, due to the mixing term {Sp\F\Ph). Thus, 

'^physical — C\ Y ph + C2 Y Sp , (36) 

with ci and C2 some constants. In other words, a spurious component Y sp is also generated which will be automatically 
projected out when computing the spectroscopic amplitude x (|2^). 

It is important to recognize that such spurious contributions are indeed needed since they account for the differences 
of the three Faddeev components (EG) . The relation between the usual Faddeev components for a given physical or 



spurious state can be inferred from the basis sets (CI). There it is shown that all the Faddeev components x*^*^ of a 
state Yp/i or Y sp are equal up, to a sign. As a consequence, if a general solution is a pure physical state Yp/j, all 
its Faddeev components cannot differ from each other in a significant way. Having a mixing between physical and 
spurious states allows the possibility of obtaining two independent Faddeev components. This result corresponds to 
the physical ingredients which involve identical ph phonons for the components x'^' and x^'^-' but a pp phonon for 

x(3). 

When all the Faddeev components are summed to generate the full x in (|2^), the contribution of the spurious states 
cancels out. Thus, for any nonspurious solution of the Faddeev equations, only the contribution from physical states 
Yp^ is needed to determine the 2plh propagator. By looking at Eq. (p5|), it is easy to see that these contributions 
can be directly obtained by diagonalizing the upper-left block 

Yj3\ = {Ph\e\Ph) Y^^ , (37) 

where to is used to label the solutions. The solutions of Eq. (pTl) are sufficient to determine the 2plh propagator. For 
some applications one may need the individual components x^^^. In that case, the contribution from spurious states 
Y Sp can be determined by solving the remaining part of the Faddeev equations 

ujm = {Sp\F\Ph) Y^\ + {Sp\F\Sp) Y^s^ ■ (38) 

We note that if the upper-right block of Eq. ( p5| ) is not zero, a mixing between the Yph and Ysp states occurs for all 
the eigenstates of the Faddeev hamiltonian. In this situation, the spurious eigenvalues will differ from the unperturbed 
energies and all of the solutions of the Faddeev equations will contain a component Yp/j. The Faddeev formalism 
would therefore become useless, since it would no longer be possible to discern between "good" and "bad" solutions. 
In Appendix ^ we show how the correct behavior of spurious solutions is related to the presence of backward-going 
contributions of the DRPA F-matrices (see Fig. |^). In case these diagrams are neglected, the spurious states Ysp 
no longer diagonalize the Faddeev hamiltonian. Such diagrams may give a small contribution to the description of 
low-lying states but they are essential to make the whole formalism presented here meaningful. As a general rule, 
when deriving expansions based on the Faddeev equations, it should be kept in mind that not all possible sets of 
diagrams can be effectively summed to all orders. Instead, one must first check the consistency of the set of diagrams 
with respect to the behavior of spurious solutions. 



V. SUMMARY AND CONCLUSIONS 



The present theoretical description of the distribution of spectroscopic strength at low energies lacks important 
ingredients for a successful comparison with experimental data. One of these ingredients is a proper description of the 
coupling of sp motion to low- lying collective modes that are present in the system. Recent calculations for ^^O [ p^ , 
for example, only include a TDA description of these collective modes. A new method is proposed here to study the 
influence of pp and ph RPA correlations on the sp propagator for a system with a finite number of fermions. This 
method is formulated in the context of SCGF theory by evaluating the nucleon self-energy in terms of the 2plh and 
2hlp propagators. The description of the 2plh (or 2hlp) excitations has been studied by using the Faddeev formalism, 
which is usually applied to solve the three-body problem. The Faddeev formalism is necessary since we consider the 
collective pp and ph RPA phonons as the basic building blocks to describe the 2plh motion. 

The computational scheme presented here employs only two-time propagators, thus leading to a tractable set of 
equations. At the same time the contributions of pp and ph RPA phonons have been consistently summed to all orders 
thereby including the physical effects that appear to be relevant for the study of the ^^O nucleus. Unlike previous 
calculations in which ph phonons have been included, the present formalism takes the Pauli-exchange correlations 
properly into account up to the 2plh level. 
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In deriving the set of Faddeev equations a formulation has been chosen which involves only a single diagonalization 
for the 2plh fragments. The appearance of spurious solutions has also been discussed in some detail, showing that the 
inclusion of the contribution of certain diagrams is necessary to separate such spurious solutions from the physically 
meaningful ones. When this separation occurs, it is straightforward to project out the physical eigenstates from the 
Faddeev equations, thereby eliminating the spurious ones. 

The Faddeev formalism has been used to include specific correlations corresponding to pp and ph phonons in a 
natural way. Extensions to the inclusion of more complicated excitations like the extended DRPA can be obtained 
in a convenient way by starting from the formalism presented in Sec. H. 

The formalism presented here appears practical for describing the spectroscopic strength in ^^O in a similar space 
as was employed in []l0| . This implementation is currently in progress and will be reported elsewhere. 
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APPENDIX A: INTERACTION BOXES. 



The r*^^'^'' matrix (14a) obtained by solving the DRPA equation has the following Lehmann representation 



(A"+ 



A"+ 



n+ 



uj — e. 



r+ 



Ar> Jlj) 



-fc- 



IT] 



(Al) 



in w hich n+ (fc— ) label the forward-going (backward-going) contributions. In obtaining the F^^^ vertex given by 
Eq. (|l3b|) we want to keep both the forward- and backward-going terms of ( Al ) . This implies that all three diagrams of 
Fig. H are included. The main problem encountered when working with dressed propagators is that the contribution of 
these three diagrams do not factorize in an expression of the form G^^FG*'^ when G°^ is represented by a propagator 
of the form 



G' 



o> 



ni ,n2 ,A; 



{el 



1 + 



(A2) 



This fa ctor ization cannot be made because of the implicit sums over the particle and hole excitation indices (n, n' , fc) 
in Eq. (A2). For example, the hole label fc cannot change in Fig. |3| [p"2|j3C| , |35|] . 



This difhculty can be overcome with a slight reformulation of the problem. We no longer regard r(3) and G°> only 
as functions of the model space indices (a, 7), but instead assign an additional dependence on the particle and hole 
indices. Thus promoting the {n, n' , fc) quantum numbers to external indices, the Lehmann representation of G°^ ( [T^ ) 
will contain at most one pole for every matrix element. As a consequence, all the components ( p^ appearing in the 
Faddeev equations have to be reformulated in the same way. The original propagators can then be retrieved at the 
end by summing the solutions over all the particle and hole fragments. With this procedure it becomes possible to 
write the sum of the three diagrams in Fig. ||in terms a matrix product of two G°^ (|l^) propagators and the following 
vertex 



F 



(3) 

fi7i^ vn^Xkx ,anal3nj3'yk-y 



(^) 



1 ^kx.k-,. 



y. 



E 

fc- 



i^k- £n„)(£fc_ ^"3) 



r- 



(A3) 
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which corresponds to the expression for the pp interaction box (13b). With this prescription, we are able to write an 
expansion that sums diagrams hke those of Fig ^. This is achieved at the cost of an increased size of the matrices 
to be dealt with. After further manipulation, it is possibl e to a void this complication by dropping the dependence 
on the model space indices (a, /?, 7), as explained in Sec. IIIB . The advantage of the present procedure lies in the 
possibility to diagonalize the Faddee v amp litudes in one step instead of solving the equations with energy-dependent 
vertex functions as discussed in Sec. Ill C . The expressions for the ph interaction boxes F*^^^ and F*^^^ are derived in 
a completely analogous way. 



APPENDIX B: DRESSED RPA EQUATIONS. 



To clarify the notations used in the paper, we give here a brief overview of the DRPA equation for the pp interaction 
matrix. We also give the explicit expressions for the normalization and closure relations used in the development of 
the formalism. 

The pp-DRPA equation is derived from Eq. (14a) by choosing K^^^^^ = Vap^-yS and is shown in Fig. |^. Using 

the Lehmann representation (Al) and extracting the poles e^^ (Sfe-) from the DRPA equation, we get the usual 
eigenvalue problem 



A 



n+(k~) 
af3 



_ .n+{k-) (0) / ,N 1 T/ 



-71 + 



(Bl) 



where n+ (fc— ) refer to the forward-going (backward-going) solutions. It is useful to introduce the following notation, 
in analogy to the convention introduced in Eq. (pO|) for the Faddeev components. 



(B2a) 



Tjk — 

ni ,712 



(B2b) 



fe- 

ki ,k2 



ykiyk2 An+ 



V2 {e^X - ^fei - £k2 



(B2c) 



(B2d) 



These represents the generalization to dressed propagators of the usual RPA components (the ^2 has been inserted 

only in the pp case for convenience). In Eq. ( p^ ) the quantities X {y) and (e^) represent the spectroscopic 

pi/ \ 

amplitudes and the poles of the forward-going (backward-going) part of the one-body propagator, while £j^_X(k--) 



the eigenvalues of the DRPA equation (Bl). 

The normalization condition for the DRPA solutions, given in terms of the components (|B2[), is the generalization 
of the normalization for the usual RPA pq] and can be put in matrix notation as 



" f/t Jf 




■ / 




'U H' 




■ / 






-/ 




J V 




-/ 



while the closure relations are given by 



(B3) 



' U H' 




■ / 




■[/t Jf 


1 




J V 




-/ 




i/t yt 


2 


jex _j 



(B4) 



where U, H, J and V are the matrices containing the elements of Eq. (B2). In dealing with the formalism for the 
Faddeev equations, it is also useful to introduce the following two matrices: 



12 



(B5a) 



k- 



V2 



A 



Q/3 



(B5b) 



which are trivially related to the com pon ents ( B2a ) and ( B2d ). 

The matrix elements given in Eqs. (B2) and (B5) correspond to the matrices i?^^^ U'^^\ and T*^'^^ introduced after 
Eq.(p3) for the 2plh Faddeev expansion (J, V, and W being the corresponding ones for the 2hlp expansion). 



APPENDIX C: PROPERTIES OF SPURIOUS STATES 

The set of solutions X_i ( ^l]) which satisfy the Pauli requirements can be divided in two subsets of physical Vph 
and spurious states Vsp- Orthogonal basis sets for these two spaces are given by 



Y 



Ph 



^ ^ex \ ^ y^^ g^^^ 



U U 



U 

ex 



u u'' 



e Vsp 



(Cl) 



where the u represent unit vectors which belong to the space (M). Their components are given by 



(C2) 



with n' , n" and k' fixed fragmentation indices which label all the possible u. The vectors u"^^ = /"^^ u are given by 
the exchange of the two particle indices n\ and 712. The physical states Yp/j are characterized by the fact that they 
do not produce vanishing spectroscopic amplitudes while the spurious states Ygp do. Thus Ysp represents the kernel 
of Eq. (||). 

It is clear from Eq. ( |Cl| ) that the Yp^ states span a space equivalent to the space of antisymmetric vectors x ([2l|), 
thus Vph = Va- Analogously, the Ygp states depend on both the symmetric and the antisymmetric parts of the u 
vectors, which implies Vsp = V)i (8) Vs. Therefore, the vectors ( plD form a basis for the full antisymmetric Faddeev 
space Vf (||). 

In general, the physical and spurious states (Cl) defined here are not solutions of the Faddeev equations (|2^), 
they simply define a basis over which these solutions can be expanded. Nevertheless, for both the normal three-body 
Faddeev equations and the expansion proposed in this paper, it can be seen that the spurious states Ygp (and only 
those) diagonalize the Faddeev hamiltonian. The eigenvalues correspond to the poles of the three freely propagating 



lines a; = £„, + e„„ — e^T (jl^). This feature serves as a sum rule on the solutions of the Faddeev equations and (unlike 
the case of 3-body systems) is not always satisfied when applying the formalism to particle and hole excitations. 
Instead this property depends on the diagrams included in the expansion and a proper set of diagrams needs to 
be employed in order to apply the Faddeev formalism. For the particular Faddeev expansion described here, this 
constraint is achieved by including the backward-going terms of DRPA phonons in the F*^') matrices and by using the 
closure relations (B4), which turn out to play an important role. In the following, the proof that the Ygp states of 
Eq. ( pl| ) actually represent a set of spurious solutions of the present Faddeev formalism is outlined. This also clarifies 
the relationship between the correct behavior of the spurious solutions and the backward-going DRPA diagrams. 

Consider a spurious state Y^^of the form (|C1|), with u given by Eq. (C2) and eigenvalue w. We now observe that 
the matrices [/ ( B2a ) and T ( B5a ) differ from each other only by an energy denominator. In particular, we have 



1 



D 



U 



(C3) 



were U, T and D are defined in Sec. III. If the eigenvalue is given by Cb 



e^, the equivalence of the 

left- and right-hand side holds only for the matrix elements having the same indices (n', n", k'). Indeed, only in that 
case the denominator a) — D will be equal to the one in Eq. ( B2a ). On the other hand, we see from Eq. (C2) that 
the components of u are nonzero only for the same indices. This allows the substitution of the — > in Eq. (|C3| ) by an 
equal sign when acting on the vector Y sp- Substituting Eq. ( |C^ ) into Eq. (|5|) and using the closure relations of the 
DRPA, we obtain the equation 



■Sp 



= MY 



Sp 



(C4) 
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which is valid only for the specific state Ysp, labeled by the indices {n',n", k'). The last equation is satisfied for a 
spurious states of the form (CI) but not for the corresponding physical state Yph- Thus, we have obtained a set of 
spurious solutions of the Faddeev equations which form an orthogonal basis of Vsp. 



In this proof, we note that the closure relation (B4) can be applied to derive Eq.(C4) because of the presence of 
the backward-going term HH^ in Eq. (^), which comes directly from the last diagram of Fig. |^. 
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FIG. 1. Example of an approximation for the self-energy. Although this approximation contains both ph and pp correlations 
it would generate incorrect results, due to the need of subtracting the 2nd order term to avoid double counting. 
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FIG. 2. Bethe-Salpeter equation (g) for the 2plh propagator with vertices given by (|9|). The irreducible interaction vertices 
for the pp and two ph channels are denoted by K^^ and if'', respectively. The irreducible vertex involving all three lines 
simultaneously is denoted by K^^^ . 



FIG. 3. Diagrams that are included in the definition of the vertex for the pp channel. Her e AF'' and AF^ are the forward- 
and backward-going part of the energy dependent contribution to the pp DRPA vertex (Al). The contribution of these three 
diagrams can be factorized in an expression of the form G"^ F'^''' G"^ only after having redefined the propagators G"^ and F'^-* 
to depend also on the particle and hole fragmentation indices {n,n ,k). The last diagram has a smaller effect on the physical 
solutions of the problem, although it is essential for the elimination of spurious solutions. 



FIG. 4. A combination of two diagrams of the type shown here can be used to connect the 2plh and 2hlp propagators. 
Diagrams like these are not included in the present approximation scheme. Nevertheless, their contribution appears in the 
normalization of spectroscopic amplitudes. Explicit time-ordering is implied in this diagram. 



FIG. 5. DRPA equation for the f'"'') matrix. 
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